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Abstract 

A search for colour reconnection effects in hadronic decays of W pairs is per- 
formed with the L3 detector at centre-of-mass energies between 189 and 209 GeV. 
The analysis is based on the study of the particle flow between jets associated to 
the same W boson and between two different W bosons in qqqq events. The ratio 
of particle yields in the different interjet regions is found to be sensitive to colour 
reconnection effects implemented in some hadronisation models. The data are com- 
pared to different models with and without such effects. An extreme scenario of 
colour reconnection is ruled out. 
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1 Introduction 



According to the string model of hadronisation, the particles produced in the process e'^e^ — >■ 
W"'"W~ — > hadrons originate, in the absence of colour reconnection, from the fragmentation 
of two colour singlet strings each of which is stretched between the two quarks from a W 
boson. In this case the hadrons are uniquely associated to a particular W and there is a 
direct correspondence between the jets formed by these hadrons and the primary quarks from 
the W boson decays. Energy-momentum is separately conserved for each of the W systems. 
However, it has been suggested that interactions may occur between the decay products of 
the two W bosons The main justification for this "cross-talk" is the relatively short 

distance separating the decay vertices of the W bosons produced in e"'"e~ annihilation 0.1 
fm) compared to the typical hadronic scale (1 fm), which implies a large space-time overlap of 
the two hadronising systems. 

The main consequence of these interactions, called Colour Reconnection (CR) effects, is 
a modification of the distribution in phase space of hadrons. CR effects are thought to be 
suppressed in the hard perturbative phase, but may be more important in the soft gluon 
emission regime Pj. While hard gluons, with energy greater than the W width, are radiated 
independently from different colour singlets, soft gluons could in principle be affected by the 
colour strings of both decaying W's. Such CR would affect the number of soft particles in 
specific phase space regions, especially outside the jet cores. 

The study of CR is interesting not only for probing QCD dynamics but also for determining 
a possible bias in the W mass measurement in the four-quark channel. CR could affect the 
invariant masses of jet pairs originating from W decays. Therefore the precision with which the 
W mass may be determined using the four-quark channel depends strongly on the understanding 
of CR effects. Events where only one W decays hadronically are unaffected by CR. 

Previous LEP studies of CR, performed at centre-of-mass energy i/i < 183 GeV, were based 
on charged particle multiplicity and momentum distributions 

The analysis presented in this paper uses the method suggested in Reference El based on 
energy and particle flow to probe the string topology of four-quark events to search for particular 
effects of particle depletion and enhancement. The results are based on 627 pb~^ of data 
collected with the L3 detector 7j at y^=189— 209 GeV. Comparisons with various models are 
made at detector level and the compatibility with the existence of CR effects in various models 
is investigated. 

2 Colour Reconnection Models 

Several phenomenological models have been proposed [3 El 131011113 El to describe CR effects 
in e"'"e~ W"'"W~ hadrons events. The analysis presented in this paper is performed with 
some of those CR models, which are implemented in the PYTHIA J2|; ARIADNE and 
HERWIG HH Monte Carlo (MC) programs. 

We investigate two models by Sjostrand and Khoze i2,J implemented in PYTHIA. They 
are based on rearrangement of the string configuration during the fragmentation process. They 
follow the space-time evolution of the strings and allow local reconnections if the strings overlap 
or cross, depending on the string definition (elongated bags or vortex lines). 

In the type I model (SKI) the strings are associated with colour flux tubes having a sig- 
nificant transverse extension. The reconnection occurs when these tubes overlap and only one 
reconnection is allowed, the one with the largest overlap volume. The reconnection probability 
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depends on this volume of overlap and is controlled by one free parameter, ki, which can be 
varied in the model to generate event samples with different fractions of reconnected events. 
The relation with the event reconnection probability (Preco) is given by the following formula: 

Preco = 1 - exp(-/A;i) (1) 

where / is a function of the overlap volume of the two strings, which depends on W-pair kine- 
matics varying with ^/s. The default value of ki is 0.6 (2j, which corresponds to a reconnection 
probability of about 30% at -y/s=189 GeV. This analysis is performed with three different values 
of ki. 0.6, 3 and 1000, corresponding to reconnection probabilities at i/s=189 GeV of about 
30%, 66% and nearly 100%, respectively. 

In the type II model (SKII) the strings have no lateral extent and the reconnection occurs, 
with unit probability, when they cross. The fraction of reconnected events in this model is of 
the order of 30% at yi=189 GeV. 

The CR model implemented in ARIADNE is based on reconnection of coloured dipoles 
before the string fragmentation takes place [H] . In the AR2 scheme, which is investigated here, 
reconnections are allowed if they reduce the string length. While reconnections within a W are 
allowed at all scales, those between W's are only allowed after the parton showers have evolved 
down to gluon energies less than 2 GeV. At ^/s =189 GeV they affect about 55% of the events. 

The CR scheme implemented in HERWIG is, as for the string fragmentation, a local phe- 
nomenon since the cluster fragmentation process follows the space-time development. In this 
model PO] the clusters are rearranged if their space-time extension is reduced. This rearrange- 
ment occurs with a probability equal to l/Ai'^oiour) with default value A^'coiour = 3, giving about 
23% of reconnected events. 

All probabilities discussed above are derived as fraction of events where at least one recon- 
nection occurs either within the same W or between two W's. 

3 Event Selection 

The energy measured in the electromagnetic and hadronic calorimeters and in the tracking 
chamber is used to select e"'"e~ — >■ W"'"W~ — hadrons events. The total visible energy (-Evis) 
and the energy imbalance parallel {E\\) and perpendicular {E±) to the beam direction are 
measured. The number of clusters, defined as objects obtained from a non-linear combination 
of charged tracks with a transverse momentum greater than 100 MeV and calorimetric clusters 
with a minimum energy of 100 MeV, is denoted by A^^duster- The selection criteria are: 

^vis/yi>0.7; EJE,,,< 0.2- \E\\\/ E^^, < 0.2- N,i^,tcr>^0. 

In addition the events must have 4 jets reconstructed with the Durham algorithm J3] with 
Vent = 0.01. To reduce the contamination from semileptonic W decays, events with energetic 
/i or e are rejected. Events with hard initial state radiation (ISR) are rejected as described 
in Reference |^ Additional criteria select events with nearly perfect quark-jet association, 
necessary for the study of particle and energy flow between jets. The two largest interjet angles 
are required to be between 100° and 140° and not adjacent. The two other interjet angles must 
be less than 100°. This selection guarantees similar sharing of energy between the four primary 
partons with the two strings evolving back-to-back and similar interjet regions between the two 
W's. The above cuts are optimized by studying MC W+W~ events at \/s=189 GeV using the 
KORALW CT] MC generator interfaced with the PYTHIA fragmentation model without CR. 
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Relaxing the angular criteria increases the efficiency but gives lower probability to have correct 
W-jet pairing due to the more complicated event topology. 

The number of selected events, the number of expected events, the selection efficiency and 
the percentage of correct pairing are given in Table H After applying all the cuts the full 
sample contains 666 events with an average efficiency of 12% and a purity of about 85% for 
e"'"e~ — > W"'"W~ hadrons. The average probability to have the correct pairing between the 
W bosons and their associated jets is estimated to be 91%. 

The background is composed of qq(7) events and Z-pair production events, in similar 
amounts. Background from semileptonic W pair decays is found to be negligible (less than 
0.3%). The qq(7) process is modeled with the KK2F MC program interfaced with JET- 
SET [in] routines to describe the QCD processes, and the background from Z-pair production 
is simulated with PYTHIA. For CR studies W-pair events are simulated with PYTHIA. All 
MC samples are passed through a realistic detector simulation jSO] which takes into account 
time dependent detector effects and inefficiencies. 

4 Particle- and Energy-Flow Distributions 

The algorithm to build the particle- and energy-flow distributions [H] (Figure starts by 
defining the plane spanned by the most energetic jet (jet 1) and the closest jet making an 
angle with jet 1 greater than 100° which is most likely associated to the same W (jet 2). For 
each event, the momentum vector direction of each particle is then projected on to this plane. 
The particle and energy flows are measured as a function of the angle, between jet 1 and 
the projected momentum vector for the particles located between jets 1 and 2. In order to 
take into account the fact that the W-pair events are not planar a new plane is defined for 
each remaining pair of adjacent jets. In this four-plane configuration the angle is defined 
as increasing from jet 1 toward jet 2, then to the closest jet from the other W (jet 3) toward 
the remaining jet (jet 4) and back to jet 1. The angle (pj^i of a particle i having a projected 
momentum vector located between jets j and j + 1 is calculated in the plane spanned by these 
two jets. A particle i making an angle 0j with respect to jet 1 adds an entry equal to 1 in the 
particle-flow distribution and adds an entry equal to its energy, normalised to the total event 
energy, in the energy-flow distribution for the corresponding bin. 

The distributions are calculated using, for the particle definition, the clusters defined in the 
previous section. 

Figure El shows the particle- and energy-flow distributions obtained for the data and the 
MC predictions at detector level by using only the ffist plane for projecting all the particles. 
The data and MC distributions agree over the full angular range in both cases. 

In order to compare the interjet regions the angles in the planes are rescaled by the angle 
between the two closest jets. For a particle i located between jets j and j + 1 the rescaled angle 
is 

0r = J - 1 + (2) 

where (pj^i is the angle between jet j and particle i and V'ij+i is the angle between jets j and 
j + 1. With this definition the four jets have fixed rescaled angle values equal to 0, 1, 2 and 3. 

Figure Ok shows the rescaled particle-flow distribution normalised to the number of events 
after a bin-by-bin background subtraction for the data and MC predictions without CR and 
for the SKI model with /ci=1000, later referred to as SKI 100%. As expected, the latter shows 
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some depletion in the number of particles in the intra-W regions spanned by the two W bosons 
(regions A and B) and some particle enhancement in the two inter-W regions (regions C and 
D) when compared to the model without CR (no-CR). 

To improve the sensitivity to CR effects the particle flows in regions A and B are averaged 
as are the particle flows in regions C and D. The results are shown in Figures Ob andEfc where 
the angle is redefined to be in the range [0,1]. MC studies at particle level with particles having 
a momentum greater than 100 MeV show that the CR effects are consistent with the detector 
level results and have similar magnitudes. 

The ratio of the particle flow between the quarks from the same W to that between quarks 
from different W's is found to be a sensitive observable to cross-talk effects as predicted by the 
SKI model. These ratios, computed from the particle- and energy-flow distributions at detector 
level, are shown in Figure ID for the data, the PYTHIA prediction without CR, the SKI model 
with fci=3 and SKI 100%. 

The differences between the models with and without CR are larger in the middle of the 
interjet regions. Therefore, in order to quantify the CR effects the ratio R is computed in 
an interval, 0.2 < ^resc < 0.8, optimized with respect to the sensitivity to SKI 100%. The 
corresponding variables for particle and energy flow are defined as follows: 



i?N = / ■ /i^-^^ d0 / / ■ /;^+° d0 and R^ = [' fi^"" d<p/ I • d0 (3) 

J0.2 / J0.2 J0.2 I JO.2 

where, in a region i, 

= and = (4) 

A^evt d0 E d(f) 

The measured values of and R^ obtained at each centre-of-mass energy are summarised 
in Table |21 Correlations in the particle rates between the four interjet regions are taken into 
account by constructing the full covariance matrix. This results in an increase of about 20% of 
the statistical uncertainty. The values obtained with the complete data sample are: 

i?N = 0.911 ± 0.023 (stat.) 
Re = 0.719 ± 0.035 (stat.) 

An estimate of the sensitivity to the SKI 100% model, shows that i?N is 2.6 times more sensitive 
than Re- Accordingly, the following results and discussion are only based on R^. 

FigureElshows the measured R^ as a function of y/s together with PYTHIA no-CR and SKI 
model predictions. The energy dependence originating from the different pairing purities and 
jet configurations is in agreement with the model predictions. For the PYTHIA SKI predictions, 
the ratio decreases with the reconnection probability over the whole energy range with similar 
magnitude. The data indicate little or no CR. 



5 Semileptonic Decays 

To verify the quality of the MC simulation of the W — qq fragmentation process and the 
possible biases which may arise when determining the particle yields between reconstructed 
jets in the detector, the particle- and energy-flow distributions are investigated in e^e~ — *■ 
W"''W~ qqlu where I = e, fi. For this analysis events are selected with high multiplicity. 
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large missing momentum and a high energy electron or muon. The missing momentum is 
considered as a fictitious particle in order to apply the Durham jet algorithm to select 4-jet 
events with ?/cut=0.01. 

The same angular criteria on the four interjet angles as applied in the fully hadronic channel 
are used here. The purity obtained after selection is about 96% and the efficiency is about 12%. 
The number of selected semileptonic events is 315 with an expectation of 314.5 events. Particle- 
and energy-fiow distributions are built in a similar way as in the fully hadronic channel with 
the additional requirement that the charged lepton should be in jet 3 or 4. Figure IHK shows the 
corresponding particle-fiow distribution projected on to the plane of jets 1 and 2 for the data 
and the KORALW MC prediction. There is good agreement between data and MC over the 
whole distribution. Figure IHb shows the rescaled particle-fiow distribution where the structure 
of the two different W's is clearly visible. The region between jet 1 and jet 2 corresponds to 
the hadronically decaying W (Wi) and the region between jet 3 and jet 4 corresponds to the 
W decaying semileptonically (^2)- The activity in the W2 region is mainly due to low energy 
fragments from the hadronic decay of the first W. A comparison of data and MC for the particle 
fiow obtained by summing the regions Wi and W2 is shown in Figure [7^. The ratio between the 
data and the MC distributions is shown in Figure UJp- This ratio is consistent with unity over 
the whole range. This result gives additional confidence in the correctness of the modelling of 
the fragmentation process of quark pairs according to the fragmentation parameters used in 
KORALW and PYTHIA as well as the particle fiow definition and reconstruction. 

In the absence of CR effects, the activity found in regions A+B of a fully hadronic event 
should be equivalent to twice the particle activity in the regions W1+W2 of the distribution 
for a semileptonic event. Figure Ub shows the ratio of the particle flow in four-quark events 
divided by twice the particle flow in semileptonic events for the data and the predictions from 
no-CR PYTHIA MC and the SKI 100% model. The CR model shows the expected deflcit 
in the hadronic channel compared to the semileptonic one. The data are consistent with the 
no-CR scheme but the large statistical uncertainty prevents a quantitative statement based on 
this model-independent comparison. 

6 Systematic Uncertainties 

Several sources of systematic uncertainties are investigated. The flrst important test is whether 
the result depends on the deflnition of the particles. The analysis is repeated using calorimetric 
clusters only. Half the difference between the two analyses is assigned as the uncertainty due 
to this effect. This is found to be the dominant systematic uncertainty. 

The second source of systematic uncertainty is the limited knowledge of quark fragmentation 
modelling. The systematic effect in the qq(7) background is estimated by comparing results 
using the JETSET and HERWIG MC programs. The corresponding uncertainty is assigned to 
be half the difference between the two models. 

The systematic uncertainty from quark fragmentation modelling in W-pair events is esti- 
mated by comparing results using PYTHIA, HERWIG and ARIADNE MC samples without 
CR. The uncertainty is assigned as the RMS between the i?N values obtained with the three 
fragmentation models. Such comparisons between different models test also possible effects of 
different fragmentation schemes which are not taken into account when varying only fragmen- 
tation parameters within one particular model. 

Another source associated with fragmentation modelling is the effect of Bose-Einstein cor- 
relations (BEG) in hadronic W decays. This effect is estimated by repeating the analysis using 
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a MC sample with BEC only between particles originating from the same W. An uncertainty 
is assigned equal to half the difference with the default MC which includes full BEC simulation 
(BE32 option) 21^ in W pairs. The sensitivity of the i?N variable to BEC is found to be small. 

The third main source of systematic uncertainty is the background estimation. The qq(7) 
background which is subtracted corresponds mainly to QCD four-jet events for which the rate 
is not well modelled by parton shower programs. PYTHIA underestimates, by about 10%, 
the four-jet rate in the selected phase space region |221- A systematic uncertainty is estimated 
by varying the qq(7) cross section by ± 5% after correcting the corresponding background by 
+5%. This correction increases the value of -Rn by 0.004. 

A last and small systematic uncertainty is associated with Z-pair production. It is estimated 
by varying the corresponding cross section by ± 10%. This variation takes into account all 
possible uncertainties pertaining to the hadronic channel, from final state interaction effects to 
the theoretical knowledge of the hadronic cross section. 

A summary of the different contributions to the systematic uncertainty is given in Table El 

The ratio obtained by taking into account the systematic uncertainties is then: 

i?N = 0.915 ± 0.023 (stat.) ± 0.021 (syst.) 

7 Comparison with Models 

The i?N values predicted by the PYTHIA no-CR, SKI, SKII, ARIADNE no-CR, AR2, HERWIG 
no-CR and HERWIG CR models are given in Tabled 

The data disfavour extreme scenarios of CR. A comparison with ARIADNE and HERWIG 
shows that the CR schemes implemented in these two models do not modify significantly the 
interjet particle activity in the hadronic W-pair decay events. Thus it is not possible to constrain 
either of these models in the present analysis. 

The dependence of on the reconnection probability is investigated with the SKI model. 
For this, four MC samples are used: the no-CR sample and those with ki=0.6, 3 and 1000. 
In the SKI model the fraction of reconnected events is controlled by the ki parameter and the 
dependence of on ki is parametrized as -RnI^i) = Pi(l — exp(— p2^i)) + Ps where pi are 
free parameters. A fit to the data is performed. The minimum is at ki = 0.08. This 
value corresponds to about 6% reconnection probability at y^=189 GeV. Within the large 
uncertainty the result is also consistent with no CR effect. 

The upper limits on ki at 68% and 95% confidence level are derived as 1.1 and 2.1 respec- 
tively. The corresponding reconnection probabilities at = 189 GeV are 45% and 64%. The 
extreme SKI scenario, in which CR occurs in essentially all events, is disfavoured by 4.9 a. 
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Table 1: Average centre-of-mass energies, integrated luminosities (£), number of selected events 
(A^events), number of expected events (A^mc), selection efficiency (e) and percentage of correct 
jet pairing (tt) for the particle flow analysis. The combined figures are given in the last row. 



(V^) (GoV) 




Re 


188.6 


0.820 ± 0.037 


0.610 ± 0.047 


191.6 


0.929 ± 0.093 


0.822 ± 0.133 


195.5 


0.948 ± 0.059 


0.774 ± 0.077 


199.5 


1.004 ± 0.067 


0.871 ± 0.095 


201.7 


0.770 ± 0.086 


0.626 ± 0.130 


205.1 


1.033 ± 0.083 


0.756 ± 0.111 


206.6 


0.958 ± 0.068 


0.781 ± 0.096 



Table 2: Measured i?N and Re values as a function of energy with their statistical uncertainties. 
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oource 


(^Rn 


Energy flow objects 


(J.UId 


qq fragmentation 


0.009 


WW fragmentation 


0.008 


BEC 


0.003 


4-jet background rate 


0.004 


ZZ background 


0.002 


Total 


0.021 



Table 3: Contributions to the systematic uncertainties on i?N- 





Rn 


Data 


0.915 ± 0.023 ± 0.021 


PYTHIA no-CR 


0.918 ± 0.003 


SKI (A;i=0.6) 


0.896 ± 0.003 


SKI (A;i=3.0) 


0.843 ± 0.003 


SKI 100% 


0.762 ± 0.003 


SKII 


0.916 ± 0.003 


ARIADNE no-CR 


0.929 ± 0.003 


AR2 


0.919 ± 0.003 


HERWIG no-CR 


0.948 ± 0.005 


HERWIG CR 


0.946 ± 0.005 



Table 4: Measured value of i?N and model predictions. 
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Figure 1: Determination of the angle for the particle i. 




Figure 2: a) Particle- and b) energy-flow distributions at ^/s = 189 — 209 GeV for data and 
MC predictions. 
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Figure 3: a) Particle- flow distribution as a function of the rescaled angle for data and for 
PYTHIA MC predictions without CR, and with the SKI 100% model. Distributions of b) 
combined intra-W particle flow and c) combined inter-W particle flow. 
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Figure 4: Ratio of a) particle- and b) energy-flow distributions (Equation 4) in regions A+B 
to that in regions C+D. Statistical uncertainties are shown. 




Figure 5: The ratio i?N as a function of ^/s at detector level for data and PYTHIA no-CR and 
SKI model predictions. The parametrisation of the energy dependence is obtained by fitting 
a second order polynomial function to the predicted MC dependence. The parametrisation 
obtained with PYTHIA no-CR gives R^{^)/Rf^{189 GeV) = -3.07x10-^5+0.1297^5-12.56. 
The dependence obtained with the SKI model 3) leads to a 2.3% change in the average 
rescaled i?N value at 189 GeV. 
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Figure 6: Particle-flow distributions a) before and b) after angle rescaling for the semileptonic 
W decays for data and KORALW prediction. 
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Figure 7: a) Particle-flow distributions as a function of the rescaled angle for the semileptonic 
W decays for data and the KORALW prediction, b) Ratio of data and MC as a function of the 
rescaled angle, c) Ratio it! of the particle flow in hadronic events divided by twice the particle 
flow in semileptonic events. 
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